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ABSTRACT 

This  study  quantifies  the  potential  impacts  on  ship  defense  high  energy  laser  (HEL)  performance  due  to 
atmospheric  effects  in  the  marine  boundary  layer  driven  by  recent  observations  and  analysis  of  worldwide  sea 
surface  temperatures  (SSTs).  The  atmospheric  effects  are  defined  using  the  worldwide  probabilistic  climatic 
database  available  in  the  High  Energy  Laser  End-to-End  Operational  Simulation  (HELEEOS)  model,  which 
includes  an  SST  database  for  the  period  1854-1997.  A  more  recent  worldwide  sea  surface  temperature  database 
was  provided  by  the  Naval  Postgraduate  School  for  the  period  1990-2008.  Mean  differences  and  trends  be¬ 
tween  the  two  SST  databases  are  used  to  deduce  possible  climate  change  impacts  on  simulated  maritime  HEL 
engagements.  The  anticipated  effects  on  HEL  propagation  performance  is  assessed  at  an  operating  wave¬ 
length  of  1.0642  /cm  across  the  world’s  oceans  and  mapped  onto  a  1°  X  1°  grid.  The  scenario  evaluated  is 
near  surface  and  nearly  horizontal  over  a  range  of  5000  m  in  which  anticipated  clear-air  maritime  aerosols 
occur.  Summer  and  winter  scenarios  are  considered.  In  addition  to  realistic  vertical  profiles  of  molecular 
and  aerosol  absorption  and  scattering,  correlated  optical  turbulence  profiles  in  probabilistic  (percentile) 
format  are  used. 


1.  Introduction 

For  the  purpose  of  evaluating  an  expected  impact  of 
climate  change  on  directed  energy  weapon  (DEW)  sys¬ 
tem  performance,  the  Air  Force  Institute  of  Technology 
(AFIT)  Center  for  Directed  Energy  (AFIT  CDE)  employed 
several  modeling  codes  to  simulate  operating  conditions. 
One  of  these  codes,  the  High  Energy  Laser  End-to-End 
Operational  Simulation  (HELEEOS),  is  perhaps  the  first 
DEW  simulation  package  to  fully  incorporate  a  corre¬ 
lated,  probabilistic  climatological  database.  The  infusion 
of  such  realistic  atmospheric  effects  into  the  simulations 
allows  HELEEOS  to  better  assess  variability  and  uncer¬ 
tainty  in  system  performance  arising  from  spatial,  spectral, 
and  temporal  variations  in  operating  conditions  (Fiorino 
et  al.  2008;  Bartell  et  al.  2005).  The  current  study  quan¬ 
tifies  the  impacts  on  ship  defense  high  energy  laser  (HEL) 
performance  due  to  atmospheric  effects  in  the  marine 
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boundary  layer  driven  by  recent  changes  in  worldwide 
sea  surface  temperature  (SST)  variations  and  trends. 
The  baseline  atmospheric  effects  are  defined  using  the 
worldwide  probabilistic  climatic  database  available 
in  the  HELEEOS  model,  which  includes  an  SST  database 
for  the  period  1854-1997.  The  Naval  Postgraduate  School 
in  Monterey,  California,  provided  AFIT  CDE  with 
a  more  recent  worldwide  sea  surface  temperature  data¬ 
base  for  the  period  1990-2008.  Mean  differences  and 
trends  between  the  two  SST  databases  are  used  to 
deduce  possible  climate  change  impacts  on  simulated 
maritime  HEL  engagements.  The  anticipated  effects  on 
HEL  propagation  performance  are  assessed  at  an  op¬ 
erating  wavelength  of  1.0642  gm  across  the  world’s 
oceans  and  mapped  onto  a  1°  X  1°  worldwide  grid.  The 
scenario  evaluated  is  near  surface  and  nearly  horizontal 
over  a  range  of  5000  m  in  which  anticipated  clear-air 
maritime  aerosols  occur.  Summer  and  winter  scenarios 
are  considered.  In  addition  to  realistic  vertical  profiles 
of  molecular  and  aerosol  absorption  and  scattering, 
correlated  optical  turbulence  profiles  in  probabilistic 
(percentile)  format  are  used. 


jamC2482 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

AUG  2010 


2.  REPORT  TYPE 


4.  TITLE  AND  SUBTITLE 

Climate  Change:  Anticipated  Effects  on  High-Energy  Laser  Weapon 
Systems  in  Maritime  Environments 

6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School, Department  of 
Oceanography , Monterey  ,C  A  ,93943 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


3.  DATES  COVERED 

00-00-2010  to  00-00-2010 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROIECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

This  study  quantifies  the  potential  impacts  on  ship  defense  high  energy  laser  (HEL)  performance  due  to 
atmospheric  effects  in  the  marine  boundary  layer  driven  by  recent  observations  and  analysis  of  worldwide 
sea  surface  temperatures  (SSTs).  The  atmospheric  effects  are  defined  using  the  worldwide  probabilistic 
climatic  database  available  in  the  High  Energy  Laser  End-to-End  Operational  Simulation  (HELEEOS) 
model,  which  includes  an  SST  database  for  the  period  185471997.  A  more  recent  worldwide  sea  surface 
temperature  database  was  provided  by  the  Naval  Postgraduate  School  for  the  period  199072008.  Mean 
differences  and  trends  between  the  two  SST  databases  are  used  to  deduce  possible  climate  change  impacts 
on  simulated  maritime  HEL  engagements.  The  anticipated  effects  on  HEL  propagation  performance  is 
assessed  at  an  operating  wavelength  of  1.0642  mm  across  the  world?s  oceans  and  mapped  onto  a  18  3  18 
grid.  The  scenario  evaluated  is  near  surface  and  nearly  horizontal  over  a  range  of  5000  m  in  which 
anticipated  clear-air  maritime  aerosols  occur.  Summer  and  winter  scenarios  are  considered.  In  addition  to 
realistic  vertical  profiles  of  molecular  and  aerosol  absorption  and  scattering,  correlated  optical  turbulence 
profiles  in  probabilistic  (percentile)  format  are  used. 

15.  SUBIECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

14 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Journal  of  Applied  Meteorology  and  Climatology  (Proof  Only) 


JOBNAME:  JAM  00#0  2010  PAGE:  2  SESS:  8  OUTPUT:  Mon  Oct  25  12:42:53  2010  Total  No.  of  Pages:  14 
/ams/j  am/0/j  amC2482 


2  JOURNAL  OF  APPLIED  METEOROLOGY  AND  CLIMATOLOGY 


a.  Description  of  the  HELEEOS  model 

The  AFIT-developed  HELEEOS  modeling  program 
includes  a  fast-calculating,  first  principles,  worldwide 
surface-to-100  km,  atmospheric  propagation  and  char¬ 
acterization  package.  This  package  enables  the  creation 
of  profiles  of  temperature,  pressure,  water  vapor  con¬ 
tent,  optical  turbulence,  atmospheric  particulates,  and 
hydrometeors  as  they  relate  to  line-by-line  layer  trans¬ 
mission,  path,  and  background  radiance  at  wavelengths 
from  the  ultraviolet  to  radio  frequencies.  Physics-based 
cloud  and  precipitation  characterizations  are  coupled 
with  physically  probable  temperature  and  moisture  ver¬ 
tical  lapse  rates  to  create  realistic  atmospheric  boundary 
layer  effects.  HELEEOS  characterizes  maritime  aerosol 
environments  using  the  Advanced  Navy  Aerosol  Model 
(ANAM)  or  various  representations  of  maritime  partic¬ 
ulates  from  the  Global  Aerosol  Dataset  (GADS).  In  the 
lowest  50  m,  HELEEOS  defines  the  maritime  optical 
turbulence  with  the  Navy  Surface  Layer  Optical  Turbu¬ 
lence  (NSLOT)  model. 

The  HELEEOS  model  enables  the  evaluation  of  un¬ 
certainty  in  near-infrared  laser  and  microwave  propa¬ 
gation  by  incorporating  probabilistic  climatological  data 
on  the  parameters  that  drive  most  major  atmospheric 
effects.  Atmospheric  parameters  investigated,  such  as 
temperature,  pressure,  water  vapor  content,  optical  tur¬ 
bulence,  and  atmospheric  particulates,  are  put  into  verti¬ 
cal  profiles  of  data  for  highly  specific  modeling  scenarios. 
Worldwide  seasonal,  diurnal,  and  geographical  spatial- 
temporal  variations  in  these  parameters  are  organized 
into  probability  density  function  (PDF)  databases  using 
a  variety  of  recently  available  resources  to  include  the 
Extreme  and  Percentile  Environmental  Reference  Tables 
(ExPERT;  Squires  et  al.  1995),  the  Master  Database  for 
Optical  Turbulence  Research  in  Support  of  Airborne  Laser 
(Bussey  et  al.  2000),  the  Global  Aerosol  Dataset  (Koepke 
et  al.  1997),  and  Air  Force  Weather  Agency  numerical 
weather  forecasting  data.  GADS  provides  aerosol  con¬ 
stituent  number  densities  on  a  5°  X  5°  grid  worldwide. 
ExPERT  mapping  software  allows  the  HELEEOS  op¬ 
erator  to  choose  from  specific  site  or  regional  upper-air 
data  to  characterize  the  correlated  molecular  absorp¬ 
tion,  aerosol  absorption,  and  scattering  by  percentile. 
The  PDF  nature  of  the  HELEEOS  atmospheric  effects 
package  enables  realistic  probabilistic  outcome  analy¬ 
ses,  which  permit  an  estimation  of  the  confidence  in  the 
calculated  probability  of  effect.  HELEEOS  users  can 
additionally  access,  display,  and  export  the  atmospheric 
data  independent  of  an  HEL  engagement  simulation 
(Fiorino  et  al.  2006).  Molecular  scattering  is  computed 
based  on  Rayleigh  theory.  Molecular  absorption  effects 
are  computed  for  the  top  13  absorbing  species  using  line 
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strength  information  from  the  High-Resolution  Trans¬ 
mission  (HITRAN)  2004  database  (Rothman  et  al.  2005) 
in  conjunction  with  a  community  standard  molecular  ab¬ 
sorption  continuum  code.  Aerosol  and  hydrometeor  scat¬ 
tering  and  absorption  are  computed  with  the  Wiscombe 
(Wiscombe  1980)-Mie  module. 

HELEEOS  provides  coverage  over  the  world’s  ocean 
regions  based  on  the  Surface  Marine  Gridded  Climatol¬ 
ogy  database,  which  is  derived  from  the  Comprehensive 
Ocean-Atmosphere  Dataset  (COADS).  This  database, 
which  includes  virtually  all  marine  data  from  1854  through 
1997  summarized  by  1°  X  1°  latitude-longitude  boxes  by 
month,  provides  the  mean,  median,  minimum,  maximum, 
standard  deviation,  and  mode  for  14  selected  elements 
from  ship  synoptic  surface  observations.  HELEEOS  cur¬ 
rently  extracts  mean  values  for  air  temperature,  sea  tem¬ 
perature,  pressure,  relative  humidity,  and  wind  speed  for 
over-ocean  locations  from  those  14  available  elements. 
The  integration  of  the  ANAM  (Gathman  et  al.  1998)  and 
the  NSLOT  (Frederickson  et  al.  2000)  codes,  both  driven 
by  parameters  from  the  Surface  Marine  Gridded  Clima¬ 
tology  database,  provides  enhanced  representation  of 
atmospheric  effects  over  all  ocean  regions  on  the  1°  X 
1°  grid. 

A  diverse  array  of  aerosol  vertical  profiles  is  also 
available.  There  are  10  profiles  defined  using  the  Optical 
Properties  of  Aerosols  and  Clouds  (OPAC;  Hess  et  al. 
1998)  code,  3  Moderate-Resolution  Atmospheric  Trans¬ 
mission  (MODTRAN)  aerosol  profiles  (Shettle  and  Fenn 
1979),  and  the  wind  speed-driven  aerosol  mixtures  from 
ANAM.  ANAM  is  the  default  aerosol  profile  for  maritime 
sites  with  the  OPAC  and  MODTRAN  definitions  avail¬ 
able  as  options. 

b.  High  energy  laser  performance  comparisons  in 
different  environments 

In  the  current  study,  “peak  irradiance’’  is  the  metric 
used  to  quantify  system  performance.  Irradiance  is  the 
measurement  of  power  per  unit  area  (W  m~2)  delivered 
to  a  target  surface.  An  increased  irradiance  value  for  a 
given  ocean  location  indicates  improved  HEL  perfor¬ 
mance.  For  this  research,  the  actual  value  of  the  peak 
irradiance  in  a  particular  environment  is  not  as  impor¬ 
tant  as  how  that  value  compares  to  the  peak  irradiance 
computed  for  the  same  engagement  setup  in  a  different 
environment.  Typically,  this  comparison  will  be  made 
through  a  ratio  of  the  peak  irradiance  for  a  climatologi- 
cally  derived  atmosphere  to  the  peak  irradiance  for  a 
standard  set  of  atmospheric  conditions. 

In  this  paper,  two  types  of  atmospheric  models  are 
compared:  a  standard  atmosphere  profile  that  does  not 
vary,  and  atmospheric  profiles  based  on  climatology  that 
do  vary  temporally  and  spatially.  The  cornerstone  of  the 
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first  type  of  profile,  the  “standard  profile,”  is  the  U.S. 
Standard  Atmosphere,  1976.  This  widely  used  model 
cannot  be  used  as  the  only  atmospheric  characterization 
for  this  study  for  several  reasons.  Foremost,  it  lacks  a 
distinct  atmospheric  boundary  layer.  The  atmosphere 
is  assigned  a  logarithmic  pressure  profile,  and  the  tem¬ 
perature  is  the  same  everywhere  at  the  surface  with 
specified  lapse  rates  for  the  primary  layers  of  the  at¬ 
mosphere.  Moisture  data  are  not  representative  at  low 
altitudes  for  most  of  the  earth’s  surface.  Also,  the  stan¬ 
dard  atmosphere  does  not  contain  aerosol  data  or  tur¬ 
bulence  profiles.  Thus,  invoking  the  standard  atmosphere 
to  quantify  FIEL  propagation  effects  can  lead  to  physi¬ 
cally  unrealistic  and  nonspatially  varying  results.  How¬ 
ever,  the  standard  atmosphere  is  well  known,  and  it  does 
allow  for  intuitive  comparisons.  For  this  study,  standard 
profiles  of  parameters  not  included  in  the  U.S.  Standard 
Atmosphere,  1976,  are  combined  with  the  standard  at¬ 
mosphere  and  are  compared  against  the  second  type  of 
profile — the  climatological  profile.  This  second  profile 
uses  meteorological,  aerosol,  and  turbulence  models  that 
are  temporally  and  spatially  variable,  allowing  for  more 
accurate  modeling  of  the  atmosphere.  These  improved 
modeling  profiles  lead  to  more  realistic  laser  weapon  sys¬ 
tem  simulations. 

c.  Establishment  of  global  synoptic  ocean 
temperature  dataset 

To  deduce  SST  departures  from  the  long-term  means 
established  in  the  143-yr  CO  ADS  dataset,  comparisons 
are  made  to  a  recently  derived  SST  dataset  that  covers 
the  period  1990-2008.  This  dataset  was  created  at  the 
Naval  Postgraduate  School  (NPS)  from  the  Global  Tem¬ 
perature  and  Salinity  Profile  Program.  Optimal  spec¬ 
tral  decomposition  was  required  to  process  the  raw 
dataset  into  a  research-usable  gridded  format,  as  de¬ 
scribed  below. 

1)  Global  Temperature  and  Salinity 
Profile  Program 

The  Global  Temperature  and  Salinity  Profile  Program 
(GTSPP)  is  a  joint  effort  of  the  International  Oceano¬ 
graphic  Data  and  Information  Exchange  committee 
(IODE)  and  the  Joint  Commission  on  Oceanography 
and  Marine  Meteorology  (JCOMM)  of  the  World  Me¬ 
teorological  Organization  (WMO)  and  the  Intergov¬ 
ernmental  Oceanographic  Commission  (IOC).  GTSPP 
creates  a  timely  global  ocean  temperature  and  salinity 
dataset  of  known  quality  in  support  of  the  World  Cli¬ 
mate  Research  Programme  (WCRP).  The  development 
of  GTSPP  began  in  1989,  and  the  program  went  into 
operation  in  November  1990  (Sun  et  al.  2009). 


Fig.  1.  A  schematic  diagram  showing  the  data  flow  and  management 
of  the  GTSPP  (from  Sun  et  al.  2009). 


Currently,  the  GTSPP  consists  of  three  components: 

1)  the  Global  Telecommunication  System  (GTS)  pro¬ 
vided  by  the  WMO  for  the  transmission  of  oceano¬ 
graphic  messages  collected  through  various  panels  in 
the  JCOMM  program;  the  GTSPP  uses  this  service  to 
acquire  real-time  data,  and  real-time  data  processing 
services  are  provided  by  the  Integrated  Science  Data 
Management  (ISDM)  service;  2)  IODE  Data  Centers, 
where  historical  data  are  acquired  either  from  other 
National  Oceanographic  Data  Center  (NODC)  offices, 
or  from  cooperation  with  projects  such  as  the  Climate 
Variability  and  Predictability  (CLIVAR)  program,  the 
World  Ocean  Database  (WOD),  and  the  Ship  of  Op¬ 
portunity  Programme  (SOOP);  and  3)  the  Continuously 
Managed  Database  (CMD).  The  NODC  provides  data- 
processing  services  for  historical  data  and  maintenance 
of  the  CMD  (also  known  as  the  GTSPP  archive).  His¬ 
torical  data  include  both  low-resolution  data  from  the 
GTS  and  the  full-resolution  data  from  expendable  bathy¬ 
thermographs  (XBTs)  or  conductivity-temperature- 
depth  instruments  (CTDs)  from  the  ships  that  provided 
the  real-time  low-resolution  data  to  the  GTS,  or  fully 
processed  and  quality  controlled  data  from  other  orga¬ 
nizations  (Fig.  1).  [FT 

The  GTSPP  data  contain  upper-ocean  temperature 
and  salinity  data  and  come  mainly  from  profiling  floats, 
XBTs,  CTDs  and  bottles.  A  dramatic  change  occurred 
in  mid-1999  with  the  initiation  of  the  Argo  Project;  this 
marked  the  beginning  of  the  use  of  temperature,  salinity, 
and  current  reports  (TESACs)  to  record  profiles  from 
robotic  profiling  floats.  It  is  clear  that  the  majority  of  the 
data  were  from  XBTs  (which  primarily  report  profiles 
with  temperature  only)  up  until  1999,  when  the  Argo 
Project  began  to  report  temperature  and  salinity  profiles 
from  profiling  floats.  During  the  period  2007-08,  the 
number  of  bathythermographs  (BATHYs)  reported 
steadily  increased  from  24  855  in  2007  to  27  775  in  2008, 
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while  the  number  of  TESACs  increased  markedly  from 
821  321  to  1  630  360  between  2007  and  2008.  A  new 
dataset  of  6869  CTD  profiles  (as  of  December  2008) 
derived  from  marine  mammals  was  made  available  for 
the  first  time  beginning  in  July  2008.  These  data  are 
useful  because  they  provide  a  high  volume  of  data  from 
areas  between  60°  and  70°S,  a  latitude  band  in  which 
data  from  other  sources  are  very  sparse. 

2)  Optimal  spectral  decomposition 

Sparse  and  noisy  ocean  GTSPP  data  need  to  be  re¬ 
analyzed  and  mapped  to  the  grid  points  before  being 
assimilated  into  numerical  models.  Any  variable  in¬ 
cluding  temperature  can  be  decomposed  into  generalized 
Fourier  series  using  the  optimal  spectral  decomposition 
(OSD)  method.  The  three-dimensional  variable  is  then 
represented  by  linear  combination  of  the  products  of  the 
basis  functions  (or  modes)  and  corresponding  Fourier 
coefficients.  If  a  rectangular  closed  ocean  basin  is  con¬ 
sidered,  the  basis  functions  are  sinusoidal  functions.  If 
a  realistic  ocean  basin  is  considered,  the  basis  functions 
are  the  eigenfunctions  of  the  three-dimensional  Laplace 
operator  with  real  topography.  The  Fourier  coefficients 
are  determined  from  observational  data  through  solving 
a  set  of  linear  algebraic  equations.  The  major  benefit  of 
using  the  OSD  method  is  that  the  boundary  conditions  for 
the  ocean  variables  (temperature,  salinity,  velocity)  are 
always  satisfied. 

Let  (x,  z)  be  horizontal  and  vertical  coordinates  and  t 
be  time.  A  variable  c(x,  z,  t)  at  depth  Zk  is  decomposed 
using  the  generalized  Fourier  series  (Chu  et  al.  2003a, b, 
2005a,b): 


data  and  therefore  are  available  prior  to  the  data  analysis. 
The  OSD  method  has  two  important  procedures:  opti¬ 
mal  mode  truncation  and  determination  of  spectral  co¬ 
efficients  {A,n}.  After  the  two  procedures,  the  generalized 
Fourier  spectrum  (1)  is  used  to  provide  data  at  regular 
grids  in  space  and  time. 

The  optimal  mode  truncation  number  Mopl  is  defined 
as  the  critical  mode  number  with  the  set  of  spectral  co¬ 
efficients  {Am}  least  sensitive  to  observational  data  sam¬ 
pling  and  noise.  For  a  sample  size  of  P  and  a  mode 
truncation  of  M,  the  spectral  coefficients  {Am\  are  esti¬ 
mated  by  the  least  squares  difference  between  the  ob¬ 
served  and  calculated  values: 


4np 


P  r 


"U)  -  X  A  Jz,  t)¥£(x,  z) 


ra= 1 


min,  (3) 


where  the  tilde  represents  the  estimated  value  at  (x,  t). 
For  homogeneously  sampled  data  with  low  noise  and 
without  systematic  error,  the  empirical  cost  function 
/emp  should  tend  to  0  monotonically  as  M  increases  to 
infinity.  The  set  of  the  spectral  coefficients  {Am}  de¬ 
pends  on  the  mode  truncation  M.  Optimal  estimation  of 
{ A  }  is  equivalent  to  the  determination  of  Mopt  (Chu 
et  al.  2003a, b).  Generally,  Mopl  is  taken  as  30-40  for  the 
basin-scale  analysis  (Chu  et  al.  2007). 

Determination  of  the  spectral  coefficients  is  achieved 
by  solving  a  set  of  linear  algebraic  equations  of  {A  (z,  /)} 
after  mode  truncation: 


M 

c(x,  Zk,t)=  A0(zk,  i)  +  X  AJzk,  t)'Vm(x,  Zk),  xe  R(zk ), 

m= 1 

(i) 

where  M  is  the  truncated  mode  number;  T,„(x,  zk)  and 
Am(z,k ,  0  are  the  orthogonal  basis  functions  (or  called 
modes)  and  the  spectral  coefficients,  respectively;  and 
R(Zk)  is  the  area  bounded  by  the  lateral  boundary  T(z(t) 
at  Zk-  The  basis  functions  {  ,l/„!(x,  z.k)\  are  eigenfunctions 
of  the  horizontal  Laplace  operator  with  the  basin  ge¬ 
ometry  and  certain  physical  boundary  conditions.  For 
temperature  and  salinity,  the  homogeneous  Neumann 
boundary  condition  is  taken  at  the  solid  boundary  T(z) 
(i.e.,  no  heat  and  salt  fluxes): 


=-A  f  ,  n  ■  V,  •'T  L  =  0,  m  =  1,  2,...,M, 

n  m  m  m  ’  h  mn  ’ 


(2) 


Aa  =  QY,  (4) 

where  a  is  the  estimated  state  vector  (L-dimensional) 
for  the  exact  state  vector  a,  A  is  a  P  X  L  coefficient 
matrix,  Q  is  a  P  X  P  square  matrix  (P  >  L),  and  Y  is  a 
/'-dimensional  observation  vector,  consisting  of  a  signal 
Y  and  a  noise  Y', 


Y  =  Y  +  Y'.  (5) 

Due  to  the  high  level  of  noise  contained  in  the  obser¬ 
vations,  the  set  of  algebraic  equations  is  ill-posed  and 
needs  to  be  solved  by  a  regularization  method.  Non¬ 
singular  orthogonal  transformation  is  conducted  through 
multiplication  of  (4)  by  a  plane  rotation  matrix  S  from 
the  left  (Chu  et  al.  2004), 

SAa  =  SQY,  (6) 


where  V);  =  d2ldx2  +  d2ldy2  and  n  is  the  unit  vector  normal 
to  T(z).  The  basis  functions  {  'I',,,}  are  independent  of  the 


which  changes  the  coefficient  matrix  and  the  source 
term  from  (A,  QY)  to  (SA,  SQY)  and  provides  the 


Journal  of  Applied  Meteorology  and  Climatology  (Proof  Only) 


JOBNAME:  JAM  00#0  2010  PAGE:  5  SESS:  8  OUTPUT:  Mon  Oct  25  12:42:56  2010  Total  No.  of  Pages:  14 
/ams/j  am/0/j  amC2482 


Month  2010 


FI  O R IN O  ET  AL. 


5 


Fig.  2.  Time  series  of  yearly  ocean  Hc  (1022  J)  for  the  upper  layer 
(0-700  m). 


opportunity  to  minimize  the  imperfection  of  the  new 
system  (6),  which  can  be  solved  by  the  usual  algebraic 
methods  such  as  the  Gauss  method. 

3)  Global  gridded  ocean  temperature 

DATASET 


More  than  6  million  temperature  profiles  (XBT,  CTD, 
Argo,  . . .)  from  GTSPP  were  transferred  from  NOAA/ 
NODC  to  NPS.  Each  profile  contains  several  hundred 
data  points.  For  any  given  month  the  original  data  size  is 
822  MB.  Our  main  objective  is  to  provide  the  scientific 
and  operational  communities  with  three-dimensional 
ocean  fields  (such  as  temperature  and  salinity)  that  have 
sufficient  resolutions  in  space  (1°  X  1°)  and  time  (1  month) 
that  can  effectively  represent  the  rapid  climate  change. 
Such  a  three-dimensional  global  synoptic  temperature 
and  salinity  dataset  was  established  from  January  1990  to 
December  2008.  For  example,  yearly  ocean  heat  content 
Hc  for  the  upper  layer  (0-700  in)  was  computed  from  the 
dataset  by 


Hc(t) 


f° 

-H  . 

x  * 

pcpT(x,  y,  z,  t)  dx  dy  dz. 


(7) 


where  (x,  y)  and  z  are  horizontal  and  vertical  coordinates, 
H  =  700  m,  and  T  is  the  yearly  temperature.  The  hori¬ 
zontal  integration  is  over  the  whole  ocean  domain.  The 
upper-ocean  Hc  has  an  evident  upward  trend  with  a  rate 
M  nearly  1.34  X  1022  J  yr  1  from  1990  to  2008  (Fig.  2). 


2.  Methodology 

Using  the  HELEEOS  model,  low-altitude  laser  en¬ 
gagements  are  simulated  to  study  atmospheric  effects 
on  the  employment  and  operation  of  military  weapon 


systems.  The  focus  in  the  current  study  is  on  SSTs  and 
marine  optical  turbulence  effects. 

HELEEOS  currently  relies  on  the  Surface  Marine 
Gridded  Climatology  database,  which  is  derived  from 
the  143-yr  COADS  dataset.  This  study  compares  the 
expected  HEL  performance  based  on  these  Surface 
Marine  Gridded  Climatology  data  and  the  new  SST  data 
obtained  from  the  NPS.  The  new  data  are  organized  by 
latitude,  longitude,  month,  and  year  and  cover  1990- 
2008.  Monthly  averages  are  calculated  over  the  19  yr  to 
generate  new  SST  means  (Fig.  3).  These  means  are  then  [F3] 
used  by  HELEEOS  to  model  the  expected  maritime 
HEL  performance. 

Engagement  and  geometry  parameters  for  the  current 
study  include  the  following: 

•  20-m  platform  altitude,  1 0  m  s  1  horizontal  velocity, 
heading  west  (270°  true); 

•  90°  relative  azimuth  to  the  target  (laser  pointing  north); 

•  5-m  target  altitude,  310  m  s  1  horizontal  velocity,  head¬ 
ing  south  (180°  true); 

•  5000-m  slant  range; 

•  1.0642-/xm  wavelength; 

•  summer  and  winter  (characterized  by  January  and 
July  data); 

•  a  Bufton  wind  profile,  which  assumes  a  10  m  s~ 1  ground 
speed  from  the  west  and  an  increase  in  speed  slightly 
with  altitude — this  defines  the  “natural  wind”;  and 

•  clear  conditions  with  no  clouds  or  rain. 

This  particular  low-altitude  engagement  scenario,  with 
its  combination  of  the  natural  crosswind  and  the  mo¬ 
tion  of  the  platform  both  perpendicular  to  the  laser  line 
of  sight  and  into  the  natural  wind,  minimizes  thermal 
blooming — essentially  “cooling”  the  laser  beam — to  high¬ 
light  the  impacts  of  marine  optical  turbulence.  Thermal 
blooming  occurs  when  energy  at  the  wavelength  the 
laser  is  operating  at  is  absorbed  by  atmospheric  con¬ 
stituents.  This  results  in  heating  of  the  beam  path  and 
a  corresponding  reduction  in  the  value  of  the  refrac¬ 
tive  index  of  the  air.  This  in  turn  results  in  a  nega¬ 
tive  lens  effect  along  the  path  and  causes  the  beam 
to  diverge.  Thermal  blooming  is  a  nonlinear  effect  and 
can  dramatically  decrease  the  effectiveness  of  a  laser 
weapon  system.  The  climatological  conditions  are  fur¬ 
ther  defined  by  the  NSLOT  turbulence  model,  heavily 
driven  by  the  air-sea  temperature  difference,  and  the 
AN  AM  model.  For  the  most  part,  the  datasets  that  are 
compared  in  this  study  are  very  similar  except  that  the 
baseline  SST  means  in  the  Surface  Marine  Gridded 
Climatology  database  (“old”  data)  are  replaced  with 
the  new  19-yr  SST  means  calculated  from  the  NPS 
data  (“new”  data).  All  other  conditions  such  as 
air  temperature  and  relative  humidity  remained  as 
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Fig.  3.  Monthly  mean  SST  calculated  from  the  gridded  GTSPP  data  from  1990  to  2008  during 
(top)  January  and  (bottom)  July. 


defined  in  the  Surface  Marine  Gridded  Climatology 
database. 

Turbulence  at  all  scales  creates  temperature  gradients 
that  cause  variations  in  the  index  of  refraction.  For 
smaller-scale  turbulence  (<100-m  length  scale,  such  as 
that  created  in  the  maritime  surface  layer  by  air-sea 
temperature  differences  and  surface  ocean  waves),  vis¬ 
ible  and  infrared  (IR)  light  will  have  different  “optical 
distances”  to  travel  over  fairly  short  paths.  This  leads  to 
changes  in  the  phase  of  the  light.  Phase  aberrations, 
when  propagating  through  space,  lead  to  changes  in  in¬ 
tensity  of  the  light  beam.  Small-scale  turbulence  as  it 
applies  to  visible  and  near-IR  light  or  laser  propagation 
is  typically  called  optical  turbulence  (Perram  et  al. 
2010).  Optical  turbulence  is  often  quantified  with  the 
index  of  refraction  structure  constant  C„.  This  param¬ 
eter  has  a  fairly  significant  wavelength  dependence 
(Fiorino  et  al.  2009);  for  example,  a  microwave  radar  at 


~10-cm  wavelength  would  measure  a  different  value 
than  a  scintillometer  operating  at  ~1  gm  would  mea¬ 
sure  for  the  same  optical  path.  This  is  due  to  the  mi¬ 
crowave  energy  being  sensitive  to  both  temperature  and 
humidity  gradients,  while  the  near-IR  light  is  mainly 
only  affected  by  temperature  gradients.  Additionally, 
focused,  near-IR  laser  energy  is  most  affected  by  optical 
turbulence  on  the  same  scale  or  smaller  than  the  laser 
aperture  diameter  (typically  less  than  1  m).  Thus,  tur¬ 
bulent  disturbances  larger  than  1  m,  such  as  those  pro¬ 
duced  by  sea  surface  waves  (Hristov  2008),  are  not  a  IAU2I 
significant  component  of  the  C'j,  parameter  for  the  near- 
IR  laser  propagation  modeled  in  this  study. 

Peak  irradiance  delivered  to  the  target  is  calculated 
for  the  old  and  new  climatological  conditions  as  defined 
above.  For  the  comparison  ratios,  a  standard  irradiance 
value  is  determined  by  HELEEOS  by  using  the  above  en¬ 
gagement  geometry  and  the  following  standard  conditions: 


Fig(s).  3  live  4/C 
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Fig.  4.  Comparison  of  (left)  optical  turbulence  strength  and  (right)  atmospheric  extinction  for  1.0642  /rm  for 
climatological  (dashed)  and  standard  conditions  (solid),  between  the  surface  and  1000-m  altitude  for  50th  percentile 
RH,  daily  average,  molecular,  and  aerosol  effects  only,  for  a  maritime  location  off  the  eastern  coast  of  Africa  (11.5°N, 
58.5°E).  The  climatological  conditions  are  based  on  NSLOT  for  optical  turbulence  and  COADS  with  ANAM 
aerosols  for  extinction.  The  standard  conditions  are  based  on  Hufnagel-Valley  5A  for  turbulence  and  U.S.  Standard 
Atmosphere,  1976,  with  MODTRAN  rural  aerosols  for  extinction. 


U.S.  Standard  Atmosphere,  1976 ,  a  rural  MODTRAN 
aerosol  model,  and  the  Hufnagel-Valley  5/7  optical  turbu¬ 
lence  model.  Hufnagel-Valley  V?  is  a  long  standing  and 
widely  used  vertical  profile  model  of  optical  turbulence, 
specifically  at  a  wavelength  of  0.55  gm,  for  which  the  Fried 
coherence  length  for  a  vertical  path  from  the  earth’s  sur¬ 
face  to  space  is  5  cm  and  the  isoplanatic  angle  is  7  gtrad 
(Hufnagel  1985).  The  Hufnagel-Valley  5/i  profile  is  de¬ 
fined  by  the  following  equation,  which  is  based  on  exten¬ 
sive  overland  thermosonde  measurements: 

C2n(h)  =  5.94  X  10“53Q^)  hwe~mm 

+  2.7  X  i(r16e-;!/1500  +  Ae~h/W0,  (8) 

where  in  this  case  h  is  the  height  (m),  W  is  the  high- 
altitude  wind  speed,  and  A  is  the  assumed  C„  surface 
value.  Typical  values  for  W  and  A,  respectively,  are 
21  m  s_1  and  1.7  X  10“ 14  nV2/3. 

The  impacts  of  the  atmospheric  boundary  layer  effects 
[F4l  on  the  results  at  1.0642  gm  are  significant.  Figure  4 
compares  low-altitude  turbulence  strength  and  extinc¬ 
tion  at  1.0642  gm  for  both  the  old  climatological  and 
standard  atmosphere  assumptions  for  a  location  off  the 
eastern  coast  of  Africa.  These  plots  provide  insight  into 
why  results  at  1.0642  gm  are  significantly  better  when 
climatological  conditions  (old  and  new)  are  assumed  for 
a  large  number  of  maritime  locations.  Climatologically 
based  turbulence  (old  and  new),  calculated  using  NSLOT 
and  driven  primarily  by  often  near-zero  air-sea  temper¬ 
ature  difference  is  two  orders  of  magnitude  less  than  the 
Hufnagel-Valley  5/i  result  at  these  low  altitudes  and 
dominates  the  results  at  the  shorter  laser  wavelengths 


for  these  extremely  low-altitude  geometries.  NSLOT  is 
a  well-verified  model  (Frederickson  et  al.  2000)  and  is 
used  in  the  simulations  to  define  the  C2  values  within  the 
first  50  m  of  altitude.  Above  50  m,  the  NSLOT  C2  value 
is  held  constant  at  the  50-m  value  until  intercepting  the 
Hufnagel-Valley  curve.  The  Hufnagel-Valley  5/i  profile 
was  developed  from  overland  optical  turbulence  mea¬ 
surements  and  includes  a  “spike”  in  near-surface  C„ 
values  to  presumably  represent  daytime  solar  heating 
effects  (Hufnagel  1985).  The  comparisons  to  NSLOT 
indicate  it  is  not  representative  of  most  near-surface 
maritime  conditions;  however,  there  is  no  accepted  mari¬ 
time  standard  optical  turbulence  model,  so  Hufnagel- 
Valley  s/7  is  used  despite  the  misfit. 

Climatologically  based  extinction  is  generally  greater 
than  that  for  the  standard  case.  At  the  shorter  wave¬ 
lengths  (such  as  1.0642  gun),  this  extinction  is  primarily 
due  to  aerosol  scattering.  The  right  panel  of  Fig.  4  dem¬ 
onstrates  how  extinction  varies  in  the  boundary  layer  in 
the  climatological  atmospheres  for  1.0642  gm.  Because 
maritime  aerosols  are  salt  based  and  hygroscopic,  they 
increase  in  size  and  cause  more  scattering  with  height 
within  the  boundary  layer  due  to  increasing  relative  hu¬ 
midity  with  height  (Fiorino  et  al.  2007).  The  decrease  of 
extinction  with  height  in  the  first  20  m  above  the  surface 
seen  in  the  right  panel  of  Fig.  4  is  due  to  wind  speed  ef¬ 
fects  on  the  multimodal  aerosol  distributions  within 
ANAM. 

HEL  performance  is  evaluated  via  analysis  of  the 
peak  irradiance  calculated  from  the  old  and  new  cli¬ 
matological  conditions  to  the  standard  irradiance  value. 
There  are  about  35  000  maritime  data  points  used  in  this 
study.  This  number  of  data  points  is  achieved  because 
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Fig.  5.  Scatterplots  of  two  locations  showing  new  SST  vs  year  with  trend  line.  The  trend  line  is  extrapolated  to 
predict  SST  for  the  years  2020  and  2050.  Uncertainty  bars  are  defined  by  the  standard  deviation  about  the  19-yr  mean 
SST  for  that  grid  point. 


both  the  Surface  Marine  Gridded  Climatology  database 
and  the  NPS  data  are  on  a  1°  X  1°  latitude-longitude 
grid. 

To  enable  a  trend  analysis  and  projection  of  potential 
SSTs,  linear  regressions  are  calculated  using  the  new 
SST  dataset.  A  simple  regression  equation  in  the  form 
of  y  =  mx  +  b  is  used  for  each  point  of  latitude  and 


longitude.  There  are  19  data  points  that  go  into  each 
calculation.  For  example,  consider  the  data  for  the  two 
grid  points  used  in  Fig.  5.  For  the  month  of  January,  IFsl 
there  are  19  temperature  measurements  at  each  of  these 
two  points  (1990-2008).  Similarly,  a  separate  regression 
is  calculated  for  the  data  in  July.  Based  on  these  regression 
equations,  SST  projections  are  made  for  2020  and  2050. 
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Irradiance  Ratio:  July.  COADS  SST  Data 
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Fig.  6.  Performance  charts  for  January  and  July  based  on  the  old  Surface  Marine  Gridded  Climatology  and  the  new  SST  data  from  NPS. 
Green  indicates  a  high-valued  ratio  that  implies  much  better  performance  in  the  climatological  conditions  than  in  the  standard  atmosphere 
conditions.  Red  indicates  a  lower-valued  ratio  that  implies  only  slightly  better  performance  in  the  climatological  conditions  than  in  the 
standard  atmosphere  conditions.  These  charts  are  based  on  the  engagement  geometry  described  above  in  section  2. 
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Irradiance  Ratio:  January.  COADS  SST  Data  Irradiance  Ratio:  July.  COADS  SST  Data 
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Fig.  7.  Histograms  of  irradiance  ratios  for  (left)  January  and  (right)  July  of  the  (top)  old  Surface  Marine  Gridded 
Climatology  calculations  and  (bottom)  new  SST  data.  Note  the  shift  to  lower-valued  ratios  in  the  histograms  on  the 
new  dataset  as  compared  to  the  old  dataset. 


The  /(-squared  values  of  the  linear  regressions  are  low 
(0.2-0.65),  but  a  linear  regression  introduces  the  least 
amount  of  bias  given  the  rather  limited  scatterplots  of  data. 
The  example  scatterplots  and  linear  regression  fits  at  two 
grid  points  in  Fig.  5  are  shown  with  uncertainty  bars  that 
are  defined  by  the  standard  deviation  about  the  19-yr  mean 
SST  for  that  grid  point.  The  projected  data  are  used  in 
HELEEOS  with  the  same  geometry/engagement  param¬ 
eters  to  model  future  (2020  and  2050)  HEL  performance. 

3.  Results 

Comparisons  are  made  between  the  irradiance  ratios 
(climatological  conditions/standard  value)  calculated 
from  the  old  Surface  Marine  Gridded  Climatology 
database  and  the  new  SST  data  from  NPS.  Over  all 
maritime  grid  points,  the  differences  are  minor,  but  the 
number  of  lower-valued  ratios  is  greater  in  the  newer 
data  from  NPS  as  compared  to  the  older  Surface  Marine 
Gridded  Climatology  data.  This  means  the  expected 


overall  HEL  performance  does  decrease  slightly  in  the 
new  SST  scenario.  Also,  the  irradiance  calculated  with 
a  standard  atmosphere  is  very  conservative  as  compared 
to  the  irradiance  calculated  with  a  temporally  and  spa¬ 
tially  varied  climatological  atmosphere  (old  or  new). 
As  a  previous  study  has  found  (Fiorino  et  al.  2008), 
the  assumption  of  standard  atmosphere  conditions — 
U.S.  Standard  Atmosphere,  1976,  with  23-km  visibility 
MODTRAN  rural  aerosols  and  Hufnagel-Valley  5h 
optical  turbulence — in  modeling  analyses  can  lead  to 
significant,  misleading  deviations  from  conclusions  de¬ 
rived  using  realistic  (climatological)  environmental 
conditions.  The  most  notable  “standard  conditions” 
deviations  from  climatology  in  this  study  are  the  lack  of 
boundary  layer  aerosol  size  and  extinction  variations 
and  the  inappropriateness  of  the  Hufnagel-Valley  5/i  C„ 
profile  in  near-surface  maritime  environments.  These 
deviations  are  apparent  in  our  data;  some  irradiance 
ratios  are  higher  than  100  in  both  the  old  and  new 
datasets.  This  means  that  the  irradiance  obtained  from 
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Irradlance  Ratio:  July  2020  SST  Data 
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Irradlance  Ratio:  January  2050  SST  Data 
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Fig.  8.  Irradiance  ratios  for  (left)  January  and  (right)  July  (top)  2020  and  (bottom)  2050.  These  irradiance  ratios  were  calculated  using 
the  same  geometry  engagement  as  described  for  the  old  and  new  climatological  plots.  They  rely  on  projected  SST  data  from  linear 
regressions  like  those  shown  in  Fig.  5. 


using  a  standard  atmosphere  can  be  as  much  as  two 
orders  of  magnitude  smaller  than  the  climatologically 
based  irradiance. 

IF6I  The  global  irradiance  ratio  charts  are  shown  in  Fig.  6 
for  January  and  July  of  the  old  Surface  Marine  Gridded 
Climatology  and  the  new  data  from  NPS.  The  top  panels 
in  Fig.  6  show  the  irradiance  ratios  for  January  and 
July  based  on  data  in  the  older  Surface  Marine  Gridded 
Climatology  database.  The  bottom  panels  in  Fig.  6  show 
irradiance  ratios  for  January  and  July  calculated  using 
the  newer  SST  data  from  NPS.  Generally,  the  results 
from  the  two  databases  follow  the  same  patterns.  Per¬ 
formance  is  better  in  the  Southern  Hemisphere  in  Jan¬ 
uary  as  indicated  by  higher  ratios,  and  in  the  Northern 
Hemisphere  in  July  for  both  databases.  This  is  primarily 
because  air-sea  temperature  differences  that  cause  more 
surface  layer  turbulence  are  generally  greater  in  winter, 
and  there  is  less  aerosol  scattering  in  summer  due  to 
lower  relative  humidity.  The  green  areas  (irradiance 
ratios  >  —90)  tend  to  be  more  broken  up  or  not  as  “solid” 
in  the  charts  from  the  newer  dataset,  but  the  overall 
change  is  very  small.  The  overall  mean  irradiance  ratio 
value  calculated  using  the  old  Surface  Marine  Gridded 
Climatology  is  —56,  and  the  mean  irradiance  ratio  value 
calculated  using  the  new  SST  data  from  NPS  is  —47  for 
a  difference  of  16%.  Modeled  performance  is  better  in 
the  standard  conditions  rather  than  the  climatological 


conditions  (irradiance  ratios  <  1)  in  only  —2%  of  the 
January  cases  (both  new  and  old),  and  —1%  of  the  July 
cases  (new  and  old).  Histograms  of  the  values  plotted  in 
Fig.  6  are  shown  in  Fig.  7.  [FT] 

By  using  the  projected  SST  data,  future  HEL  perfor¬ 
mance  is  modeled  using  HELEEOS.  The  irradiance 
ratio  charts  for  2020  and  2050  are  shown  in  Fig.  8,  and  [re] 
the  difference  between  the  2020  and  2050  levels  of  per¬ 
formance  and  the  current  level  of  performance  is  shown 
in  Fig.  9.  The  2020  and  2050  performance  comparisons  to  [re] 
current  conditions  are  made  by  subtracting  the  old  cli¬ 
matological  irradiance  ratios  (top  plots  in  Fig.  6)  from 
the  2020  or  2050  irradiance  ratios.  Figure  10  shows  his-  [FTol 
tograms  of  the  2020  irradiance  ratio-old  data  irradiance 
ratio  differences  for  January  and  July.  Only  slight  dif¬ 
ferences  are  discernible  in  the  performance  for  2020  as 
compared  to  the  current  performance,  as  is  evidenced  by 
the  large  amount  yellow,  near-zero  difference  areas  de¬ 
picted  in  the  top  panels  of  Fig.  9,  and  the  near-zero  peaks 
in  the  histograms  of  Fig.  10.  However,  when  the  gridpoint 
trends  are  projected  another  30  yr  to  2050,  the  compari¬ 
son  to  current  conditions  shows  the  expected  perfor¬ 
mance  degradation  over  large  areas  of  the  North  Atlantic 
and  Bay  of  Bengal  in  July,  and  the  South  Atlantic  in 
January.  This  is  evident  in  the  red-colored  regions  seen  in 
the  bottom  panels  of  Fig.  9.  It  is  important  to  note  that 
these  performance  changes  are  based  on  the  assumption 
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Irradlance  Ratio  Difference: 
July  2050  Data  -  July  COADS  Data 
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Fig.  9.  Irradiance  ratios  differences  for  (left)  January  and  (right)  July.  These  irradiance  ratio  differences  were  calculated  by  subtracting 
the  irradiance  ratio  value  from  the  old  climatological  database  from  the  (top)  2020  or  (bottom)  2050  irradiance  ratio.  In  general,  a  yellow 
color  indicates  little  change  in  expected  performance  from  current  conditions  to  the  projected  year,  green  indicates  improved  expected 
performance,  and  red  indicates  degraded  expected  performance. 


that  the  ocean  temperatures  will  continue  to  follow  the 
same  trends  at  each  point  for  the  next  1CM-0  yr. 

Roughly  60%  of  the  irradiance  ratios  appear  as  yellow 
in  the  top  (2020)  panels  of  Fig.  9,  a  ratio  difference  of  25 


or  less.  Most  areas  show  little  discernable  change.  Some 
areas,  such  as  the  North  Pacific  in  January,  the  Korean 
Peninsula  in  July  (Fig.  11,  right),  and  some  coastal  re¬ 
gions  around  Asia  and  Africa  in  July  do  show  a  slight 
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Irradiance  Ratio  Difference: 
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Fig.  10.  Irradiance  ratio  differences:  2020  performance  —  current  (old  data)  performance  with  histograms.  The 
difference  was  found  by  subtracting  the  performance  calculated  using  the  Surface  Marine  Gridded  Climatology 
database  from  the  predicted  performance  in  2020  based  off  of  the  projected  SST  data  from  NPS.  The  histograms  of 
the  irradiance  difference  show  over  7000  data  points  are  centered  around  zero. 
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Fig.  11.  Zoomed-in  views  of  (left)  South  America  in  January  and  (right)  the  Korean  Peninsula  in  July  showing  the  difference  between 
the  irradiance  ratio  for  the  year  2020  and  now.  Note  the  poorer  than  expected  performance  off  the  coast  of  South  America  but  the  better 
than  expected  performance  around  the  Korean  Peninsula.  These  plots  are  zoomed-in  views  of  the  top  panels  in  Fig.  9. 


improvement;  however,  many  areas  show  little  change. 
The  eastern  coast  of  South  America  in  January  (see 
Fig.  11,  left)  and  the  Atlantic  and  Indian  Oceans  in 
July  show  the  largest  decreases  in  performance.  The 
irradiance  ratio  differences  in  these  regions  constitute 
most  of  the  nonzero  values  shown  in  Fig.  10. 

The  main  driving  force  behind  the  changes  in  HEL 
performance  is  the  air-sea  temperature  difference.  Air- 
sea  temperature  difference  is  the  primary  variable  that 
determines  the  surface  layer  turbulence  that  can  dra¬ 
matically  affect  laser  propagation  in  the  lower  atmo¬ 
sphere  near  the  sea  surface.  Air  temperature  data  are 
available  in  the  Surface  Marine  Gridded  Climatology 
database  and  are  used  by  HELEEOS  in  calculations. 
However,  air  temperature  is  not  varied  in  this  study. 
Therefore,  any  change  in  SST  directly  results  in  a  change 
in  the  air-sea  temperature  difference.  These  changes  in 
air-sea  temperature  difference  either  increase  or  de¬ 
crease  deleterious  atmospheric  effects  on  expected  HEL 
performance.  This  identifies  the  main  reasoning  for  not 
allowing  the  air  temperature  to  vary  in  response  to  the 
ocean  surface  temperature  changes — it  maximizes  the 
IF12I  differences  to  simulate  the  worst-case  scenario.  Figure  12 
shows  the  expected  air-sea  temperature  differences  for 
January  and  July  2020  with  this  assumption  applied. 

Average  SST  was  calculated  for  the  Surface  Marine 
Gridded  Climatology  database  and  for  the  new  data  from 
NPS.  The  global  SST  average  for  the  Surface  Marine 
Gridded  Climatology  database  was  found  to  be  17.3°C. 
The  global  SST  average  for  the  new  data  from  NPS  was 
found  to  be  16.1°C.  It  should  be  noted,  however,  that  the 
COADS-derived  Surface  Marine  Gridded  Climatology 
database  spans  a  much  longer  time  period  (143  yr)  and 
includes  more  observations  than  the  new  dataset  from 
NPS,  which  spans  only  19  yr.  To  find  these  averages, 
only  January  and  July  of  each  dataset  were  used.  In  the 
new  dataset,  there  are  a  few  points  that  are  significantly 
cooler  than  the  COADS  data  for  the  same  points,  thereby 


lowering  the  global  mean,  but  most  points  are  slightly 
warmer  and  exhibit  trend  lines  with  small  positive  slopes, 
as  shown  on  the  left  in  Fig.  5. 

4.  Conclusions 

The  overall  conclusion  of  this  study — where  antici¬ 
pated  SST  change  is  used  to  drive  marine  surface  layer 
optical  turbulence  change  that  could  impact  maritime 
high  energy  laser  engagement  scenarios — is  that  the 
effects  of  the  anticipated  SST  changes  would  in  gen¬ 
eral  be  minor  or  negligible  in  most  oceanic  locations. 
The  analysis  of  the  COADS  dataset  versus  the  19-yr 
NPS  dataset  did  indicate  a  very  slight  decline  in  ex¬ 
pected  HEL  performance  in  the  scenario  studied  at 
a  number  of  locations,  but  most  areas  showed  little 
discernible  change. 

Analysis  of  the  COADS  data  versus  the  NPS  SSTs 
does  show  that  the  air-sea  temperature  difference  has 
become  slightly  more  negative  on  a  world- wide  scale.  At 
most  points,  this  means  SSTs  are  slightly  warmer  in  the 
19-yr  database  than  in  the  143-yr  database.  Notably,  the 
change  in  SSTs  reduces  the  air-sea  temperature  differ¬ 
ence  in  some  areas  (e.g.,  around  the  Korean  Peninsula  in 
July),  thus  leading  to  a  slight  improvement  in  antici¬ 
pated  HEL  weapon  system  performance  based  on  the 
NPS  data  over  the  COADS  dataset. 

For  the  portion  of  this  study  where  the  differences 
between  the  COADS  and  NPS  datasets  are  analyzed  for 
trends  and  projected  to  2020  and  2050  scenarios,  the  vast 
majority  of  the  oceanic  locations  show  the  changes  in 
expected  HEL  performance  to  be  insignificant.  Specifi¬ 
cally,  the  projection  analysis  indicates  areas  in  the  South 
Atlantic  to  have  the  largest  decrease  in  expected  per¬ 
formance,  while  the  North  Pacific  shows  the  greatest 
increase  in  expected  performance  in  2020  and  2050. 
The  histograms  show  that  the  frequency  of  lower-valued 
irradiance  ratios  increases  for  2020,  suggesting  more 
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Fig.  12.  Air-sea  temperature  difference  for  (left)  January  and  (right)  July  2020.  (top)  The  air-sea  temperature  difference  plots  show 
that  the  majority  of  points  exhibit  little  change  (green  shades),  (bottom)  The  histograms  show  that  most  points  indicate  a  slightly  more 
negative  air-sea  temperature  difference,  meaning  that  ocean  temperature  increased  slightly  at  most  grid  points. 


locations  experience  a  slight  decrease  in  anticipated 
HEL  system  performance  than  an  increase.  On  a  global 
scale,  the  air-sea  temperature  difference  continues  to 
become  more  negative,  meaning  sea  temperatures  in¬ 
crease  slightly.  However,  air  temperature  was  not  varied 
in  this  worst-case  sensitivity  study,  so  any  change  in  SST 
will  directly  result  in  a  change  in  the  air-sea  temperature 
difference. 
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